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ABSTRACT 

D , 

^ . Deep Inelastic scattering experiments using transversely polarised targets yield informa- 



tion on the structure function g2. By means of a free-field analysis, we study the operator 
structure of g2 and demonstrate the need for retaining the twist three mass terms in order 
to maintain current-conservation. We show that the structure function gr = 91 + 92 has a 
much simpler operator structure as compared to g2, in spite of the fact that, like g2, gr has a 
twist-three component. We demonstrate factorisation of the hadronic tensor into hard and 
soft parts for the case of gx- We show that the first moment of the gluonic contribution to 
gT vanishes, and discuss possible physical applications. 
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The recent transversely polarised deep inelastic scattering (DIS) experiments [|l| have 
opened up new avenues to explore the polarised structure of the proton. With more data, 
better accuracy and new proposed experiments [0, it would soon be possible to extract the 
twist-three contribution for the first time. Though extraction of the twist-three contributions 
is in general quite difficult P, ^, in these experiments it is possible to kinematically eliminate 
the leading twist contribution p. 

The polarised structure of the proton is characterised by two structure functions gi{x, Q"^) 
and g2{x, Q"^) which can be measured in a polarised lepton-proton DIS experiment i{k)P{p) 
— s> i{k')X{px). The spin dependent part of the proton tensor is parametrised as 

W,.{x,Q^) = ^e,,x.q^\s''(9iix,Q^)+g2ix,Q^))-^p^g2ix,Q')], (1) 
p ■ q [ ^ ^ p ■ q J 

where is the spin vector of the proton and is normalised as = — with s ■ p = 0, M 
being the target mass. The spin-dependent cross-section ^ is given by 
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where y = p ■ q/p ■ k, k = 1 + Ax'^M'^/Q'^, (p is the azimuthal angle and a is angle between 
the spin vector s and the incoming lepton momentum k. In a longitudinally polarised ex- 
periment {a = 0), the dominant contribution comes from the structure function gi{x,Q'^) 
while g2{x,Q'^) is suppressed by a factor M'^/Q'^, thus enabling the extraction of gii^x^Q"^). 
The longitudinally polarised DIS process has been studied quite extensively 0] and there 
is a considerable amount of data on gi^x.Q"^) In contrast, the extraction of g2{x,Q^) 
requires transversely polarised proton {a = 90) and further this cross-section is suppressed 
by a factor M/ \fQ^ relative to the longitudinal case. Note that at the cross-section level the 
transverse asymmetry measures the twist-three contribution while the longitudinal asym- 
metry measures the twist-two contribution. Hence the extraction of g2[Xi Q^) is much more 
complicated as compared to gi{x,Q'^). Recently, experimental information on g2{x,Q'^) has 
become available , but the data have large errors and do not provide a definite answer to 
the question of the validity of the sum-rules associated with g2{x,Q^), like the Burkhardt- 
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Cottingham (BC) sum-rule the Wandzura-Wilczek (WW) sum-rule |jT^, or the recently 



proposed Efremov-Leader-Teryaev (ELT) sum-rule [|TT|. contribution etc. 

In transversely polarised DIS experiments, the asymmetry that is measured is the virtual 
photon absorption asjTiimetry 

At n2\ y/Q^ giix,Q^)+g2ix,Q^) 

^^^^'^^ = F^M^) ' 

where Fi(x,Q'^) is the spin-averaged structure function. We see from eq. |^ that the asym- 
metry is proportional not to g2 alone, but to Qt^x, Q"^) = gi{x, Q^) +g2{x, Q^). In this letter, 
we show that the quantity gx admits of a much simpler description than does 5^2 • We suggest 
that this may help in going some way towards a fuller understanding of the transverse spin 
structure of the nucleon. We begin by discussing the free field theory analysis [0 in order to 
elucidate the operator structure of the structure functions and demonstrate the importance 
of the mass term in maintaining gauge invariance of the hadronic tensor. We then discuss 
the first moment of grix, Q^) and its relation to the spin content of the proton, and study 
the gluonic contribution to the first moment, using the Factorisation Method (FM) |]T3| . 
The hadronic tensor W^y{p, q, s) has the form 

W,.{p,q,s) = l^Jd'^ e''-^ {ps\ [UO, MO)] \ps)c , (4) 

Retaining the dominant contribution in the light-cone limit, ^"^ —>■ 0, identified as the most 
singular part of the time-ordered product of these currents on the light-cone, we find 

W,.{p, q, s) = ^e^.Ap / d'^ e''^-^ ^ S^'\e) e{^o) {ps\ : O'^^, 0) : \ps), , (5) 

where 

o'Ai^, 0) = mi'i^m + mi'i^m , (6) 

To arrive at the above result we used 

iS({,o) = -{o|r(^(0^^(o))|o) , (7) 

where order m terms are neglected. The importance of these terms will be shown later. To 
find the dominant contribution coming from these operators, we still have to make them 
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local and then pick up the dominant part. Hence, 



(9) 



As we can see, the above local operator is symmetric in fii....fin but has no definite symmetry 
in the permutation of p with any of the other indices. The dominant part of the above 
operator can be obtained, in the usual twist analysis, by decomposing into symmetric and 
mixed symmetric parts. 



{p 



(10) 



Pip.2---Pn} 12+1 Pl\P2---P-n} 

The fully symmetric part is twist-two and the mixed symmetric part is twist-three. (As 
usual, { }, [ ] mean symmetrisation and antisymmetrisation respectively). 

Let us now compute the twist-two contribution to W^y{j),q,s). Expanding the operator 
matrix element in terms of the vectors available in the theory, the most general expression 
which is fully symmetric can be written as 



{p 



^J■lP2■■■^J■n} 
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n-l 
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(n + 1)! 

where Bn{p^) are unknown scalars which contain all the non-perturbative information. 

To perform the integration in eqn.(|^) using the delta function, we define the function 
giy) such that 



h{y) 



1 ^ B^{p^] 
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Substituting eqn.(|lT]) in eqn.(^ and using the above Fourier decomposition we can perform 
the integrals and the result is 

d 
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sP-ppllA [l + x- 
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(13) 



Comparing the above expression with eqn.(p, we find 
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(15) 
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where the superscript 2 denotes the twist of the operators contributing to the structure 
functions. Note that current conservation is maintained, i.e q^W^v = 0. Let us now compute 
the twist-three contribution to W^^ip, q, s). Again using symmetry arguments we find that 

= is^p" - syn ■ pr~' ■ (16) 

Substituting this in eqn.(||) and using 
we find that 



q^('jlpP^s^]^+php(l-x4- 
\p.q J ax \ ax . 



d{x) , (18) 



where the superscript on the hadronic tensor denotes the twist. From the above equation it is 
clear that the second term does not satisfy the current conservation relation. It is important 
to realise that this non-conservation does not manifest only for the hadronic matrix elements, 
but continues to hold even if we were to compute the eqn.(^) between quark states. Doing 
this, we find 

W,. = ^^^.uxps" {{p + g)\(g° + p')6{l - x) + {q - p)\{q' - p')6{l + x)) . (19) 

Notice that current conservation is violated even at this level. It can be maintained if we 
include the mass term which we dropped in the expansion of the time-ordered product. The 
mass term turns out to be 

Adding this term to the equation (|^) and using the equation of motion for the quark fields, 
we get the manifestly current conserved form: 



»7r 



l^^,uxpq'j rfV-«5(e')e(eo)(p,^| : ^"(^,0) : \p,s) . (21) 



We find that the mass term exactly cancels the current non-conserving part appearing in 
the eqn.(^ to reproduce the above current conserved equation. The above analysis shows 
us that when we work at a given order in the twist expansion, it is important to keep both 



5 



the singular terms and the regular terms that can contribute at the given order. The final 
result that we have when we combine the twist-2 and the twist-3 contributions is as follows: 



92[X) 



d , , , p f, d\ 

—X— oix) 1 + X— 

ax p.q \ ax J 



d{x) 



1 + x 



d_' 
dx 



b(x) — d(x) 

dx 



(22) 
(23) 



Here, one can easily see that gi{x, Q^) and g2{x., Q"^) are related by WW sum-rule if the higher 
twist terms such as terms proportional to and d{x) are neglected. /^From the expression 
for g2{x,Q'^), it is interesting to note that the twist-two part b{x) and the twist-three part 
d{x) contribute to the cross-section at the same order in M/ y^Cp. Since they appear at the 
same order, it is very difficult to disentangle these operators and, in general, a measurement 
of g2{x) will be sensitive to both twist-two and twist-three operators at the same level. While 
twist-two operators have a simple parton model interpretation, higher twist operators cannot 
be described with the same simple picture. As a consequence, g2{x, Q"^) is not amenable to a 
parton model interpretation. We shall see, in the following, that there is yet another problem 
with the interpretation of g2{x), viz., the usual factorisation of the hadronic tensor into hard 
and soft parts, wherein there is a cancellation of the infrared and collinear singularities, 
does not occur for the case of (?2- On the contrary, if we consider gT instead of g2, then 
factorisation can, indeed, be demonstrated. In what follows, we proceed to verify these 
statements, using the factorisation method. 




P 



Fig. 1. Born diagram. 



The factorisation theorem ensures the separation of long distance (soft) effects from 
the short distance (hard) effects and hence in the DIS limit the quark and gluonic contribu- 



6 



tion to the polarised hadron tenor W'^'' can be factorised as 

i ''^ y 

+ f -fA,/p{y. 1^') H;:%q, yp, asili')) , (24) 
Jx y 

where i runs over the quark flavours and fi is the factorisation scale which deflnes the 
separation of short distance from the long distance part. The soft effects are contained in 
the parton distribution functions /ao/p which are proton matrix elements of certain gauge 
invariant bifocal operators made out of parton fields such as quarks and gluons. The hard 
scattering coefficients (HSC), i?^*" are perturbative and the factorisation theorem in the DIS 
limit ensures that they are free of any infrared (IR) and coUinear singularities and do not 
depend on the properties of the target. This target independence can be used to advantage: 
the HSCs can be computed order by order by replacing hadron states by asymptotic parton 
states. The contribution to various structure functions can be extracted by using appropriate 
projection operators. 





Fig. 2. Photon gluon fusion diagram 
Let us begin by evaluating the HSCs to the quark sector to leading order. Replace the 
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proton by quark target (P q) and retain terms up to order 0{a^g) in eqn.(^) 

i -^^ y 

where the superscript in the above equation denotes the order of strong couphng a^. W^l^''^*'^ 
to 0{a^) gets contribution from the Born diagram 7*g —>■ q (Fig. 1). For quarks as target, 
the distribution function can be calculated from the operator definitions (see below). To 
0{a^), f^l/g oc 6{1 — z) and hence H^^J'^''^ is same as the Born diagram. This is the usual 
statement that to leading order the parton model (PM) and the factorisation method which 
is a field theoretical generalisation of the PM are equivalent. 

We now discuss factorisation at the next-to-leading order. To evaluate H'^^f in eqn.(p^), 
we replace P ^ g, i.e 



wi^y^{x,Q^) = e / -/S/,(?/,/^^X^'^*ng,m/^^«.(/.^)) 

ix y 



+ [ -/Ai/.(2/'/^') 4l^'"*n9,m/^',«,(/.2)) , (26) 



The LHS is the subprocess cross-section ^9 ~^ (Fig- 2) and the RHS has two parts viz. 
the quark and gluon sector. We have shown that Hjj^J''^*'^ is the subprocess 7*g q. Hence 
we need to evaluate /a^/^ and /^j/g from the definitions given in eqn. (|27|j28|) , by replacing 
P ^ ^. To C(a°), oc (5(1 -z). To evaluate H^y^, we have to evaluate Wj^J^^*'^ and 

^Aq/g- '^^^ above procedure works fine for the unpolarised structure functions Fi^2{x,Q'^) 
and the longitudinally polarised structure function gi{x,Q^) |T^, |1^], but the trans- 
versely polarised structure function g2{x,Q^) turns out to be an exception. Projecting the 
contribution to g2{x,Q'^) in eqn.(^), it turns out that ifg"^'^ = 0. This is expected as we 
know from the PM that g2{x, Q^) = to leading order. As a consequence it turns out from 
eqn. (PUj) that 14^2^''''*^ ^ = -^^2^"*'^ ^- Hence the usual cancellation of IR and collinear singu- 
larities between the subprocess cross-section and the appropriate parton matrix element, to 
give a HSC free of these singularities does not seem to occur in the case of g2{x,Q'^). The 
reason for this is that the above expression is incomplete as far as the extraction of g2{x, Q"^) 
is concerned. In fact the operator structure for the g2{x, Q"^) is much more complicated than 
that of gi{x, Q^). Also, the simple minded convolution may not work for g2{x, Q"^). This has 
earlier been demonstrated using free field analysis. 
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Let us now demonstrate the claim of factorisation made for the case of gr- To leading 
twist, gT{x,Q^) gets contribution from the parton distributions functions, viz. 



/Ag/p(a;,/i^ 



/Ag/p(x,/i^ 



Air 



^fiuXa s^p" I di e 



-{vsi.\F^'mtF^\n\vsL). 



(27) 



f28) 



where the hght-cone variables have been used to denote any four vector ^'^ = (^^,^~,^t), 
with = (^° ± C,^)/V2. F^^" is the gluon field strength tensor and = V exp[ig d(~ 
y4+(C^)]^ is the path ordered exponent which restores the gauge invariance of the bilocal 
operators. We do not consider the other twist-three gluonic operator ||T^ that contributes to 
DIS as they are suppressed by strong coupling. For parton targets the above distributions 
are normalised as 



f{Aq+Aq)/a{h){z) 
fAg/a{h){z) 



h 5{l - Z) 6a,g , 



(29) 
(30) 



where /i = ±1 is the helicity of the incoming parton and z is the sub-process Bjorken variable. 

Using the procedure discussed above, we evaluate the HSCs to transversely polarised 
structure function gxi^x, Q"^) and study its factorisation properties. To 0(a°) we can project 
the gxiXjQ"^) contribution from eqn.(p5|). The LHS is the Born diagram 'y*{q)q{p) — > q{p') 
(Fig. 1) and its contribution to gx^x^Q"^) ^ 0. Using the normahsation condition eqn.(^) 
we find 



(0),7*<? 



5(1 



(31) 



From eqn. (^4]) it is clear that to leading order gT{x,Q ) gets contribution from the parton 
distribution eqn. (|27|) . At next to leading order gT{x,Q'^) gets contribution from the other 
parton distribution eqn. (^8D . To evaluate the corresponding HSC H^^'^ ^(x,Q^), we use 
eqn.(PBD and project out the grix, Q"^) contribution. This involves the calculation of the ma- 
trix element eqn.(^) between gluon states /^g/g and the sub process cross-section w!^ 
both to 0{as). 
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The sub-process cross-section W^^],^''^*^ involves the '^*{q)g{k) — > (l{p)<i{p') fusion process 
(Fig 2). This diagram is free of UV divergence but has a mass singularity, which appears 
at small scattering angles in the massless limit. We could regulate this by keeping either 
the quark or the gluon mass non-zero, but we choose to keep both particles massive as 
the prescription dependence would be explicit in this case. Projecting the contribution to 
gT{x,Q'^), by using the appropriate projection operator, we get 



(1),7*9 



2 '^s 



dL 



An J-i (1 - L'^uj-^k)^ 



-1 + z + —z 

1^ ) 



m 



1 - L^ujy - 4— z(l + L^uj^) 



, a'^ 2 r 2 2 



-2{l ~ z) + a'^^^^z - {I - z){l 



LV) - -Ml 



2, ,2n 



where tu^ = 1 — Aw? / s, s = {q + kY and L = cos9, 9 being the centre-of-mass scattering 
angle. Performing the two body phase-space integral in the centre-of-mass frame, we get 

2 as k'^z{l — zY 



(1),7*9 



277 m? — k'^ziX — 2) ' ^ ^ 

Note that the contribution to gT{x,Q'^) is independent of the InQ^ term. Both gi{x,Q^) 
and g2{x, Q"^) separately depend on the InQ^, but the combination gri^x, Q'^) is independent 
of the InQ^ term. The constant piece depends on the choice of regulator and hence is 
prescription dependent as is clearly seen. 

7±75 7±75 





Fig. 3. The 0{as) contribution to the matrix fAq/g- 

The matrix element f^l/g (Fig- 3) is evaluated using the parton distribution eqn.(p^) 
in the light-cone gauge with the replacement P g. We keep the quarks and gluon off 
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mass-shell as we did for the evaluation of cross-section. Noting that this matrix element 
is superficially divergent, we compute it using dimensional regularisation method, and the 
matrix element turns out to be 



fW -2a L 

J(Aa-\-Aa)/o " 



As is clear, the integral is convergent and reduces to a simple form 



+ k'^z{l - z) + J:^P\ 



(33) 



.(1) ^ a. k''z{l - zf 

/(Ag+A,-)/g m2 - A;2^(l - z) ' ^ ' 

The term {d — 4:)/{d — 2) in the integral gives non- vanishing finite contribution. We also 
checked the correctness of our result in the Pauli-Villars (PV) regularisation scheme. We 
could reproduce the same result in this scheme also confirming that our finite result is 
UV scheme independent. In the PV regularisation, since the integral is performed in four 
dimension, the {d — A)/{d — 2) term is absent. The analogous term comes from the integral 
with m replaced by M (PV regulator) in the limit M goes to infinity. Hence, our result is 
independent of UV scheme. This is not true in the case of operator matrix elements which one 
encounters in the evaluation of the structure functions F2{x,Q'^) and 5'i(a;,Q^) [l^. Recall 
that the matrix elements appearing in the evaluation of the QCD corrections to F2{x,Q'^) 
and (7i(x, Q^) are UV renormalisation scheme dependent. In other words, those operators 



are defined/renormalised in a definite UV renormalisation scheme say MS or momentum 
subtraction scheme or Pauli-Villars scheme. In our case, since the matrix element is finite 
the result is UV renormalisation scheme independent to this order. Observe that the masses 
we introduced to avoid IR singularities lead to two different results when one considers 
the cross-section and the matrix element separately. That is, both the cross-section and 
the matrix element are dependent on the order in which the masses go to zero. This is 
the usual prescription dependence one encounters in massless theories. The prescription- 
dependent structure of the above equation is the same as that of Substituting for 
the normalisation condition eqn. (|30D and eqn.(^,^,0) in eqn.(p6D, we get 



H^'^'"''' = . (35) 

Note that there is a cancellation of the prescription dependent pieces confirming the fac- 
torisation. Also, it turns out that the next to leading order HSC, H!j^^'"' ^ is zero and hence 
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twist-three distribution eqn.(|28|) does not contribute to griXjQ"^) to any of the moment to 
this order. The above analysis proves that the first moment of the gluon coefficient function 
is zero in FM. 

One of the important outcomes of the demonstration of factorisation for qt is that it 
admits a description in terms of a process-independent universal distribution. This distri- 
bution is no longer a parton distribution in the usual sense of the term, because of the 
twist-three contribution to Qt- But the process-independence is still useful, so that once the 
non-perturbative distribution associated with qt has been extracted in one experiment (in 
DIS, for example), it can be used to make predictions for other processes, like Drell-Yan. 

The above analysis also has interesting consequences for the first moment of the structure 
function grix, Q^). At the leading order, one would be led by the validity of the BC sum- rule 
to conclude that the first moment of griXiQ"^) is same as that of gi{x^Q'^). Though these 
first moments are measured in completely different experiments {gT{x,Q'^) in transversely 
polarised DIS and gi{x,Q'^) in longitudinally polarised DIS), they should coincide numeri- 
cally. Note that because the BC sum- rule is valid at one loop order one would expect 
that the first moment of the gluonic coefficient vanishes in the FM. Our analysis of gx using 
the factorisation method confirms this expectation. The first moments of both gi{x, Q^) and 
griXjQ"^) in FM, are related to one and the same matrix element {psltjj ^'j5ip\ps) which is 
Lorentz invariant. The argument used here is the same as the rotational invariance argument 
that may be used to justify the BC sum-rule |^. 

A related issue is that of the hard gluonic contribution to the first moment of gi via 



the anomaly [jT9|. This gluonic contribution induced through the anomaly is, in fact, a 
possible explanation for the surprisingly small value for the first moment of gi measured in 
experiments P|. In the FM JlSl, however, this contribution vanishes as long as the quark 
distributions are related to matrix elements of the standard quark field-operators which 



appear in the operator product expansion 16 1. In a parton model computation of the hard 



gluonic contribution due to the anomaly, the gluonic coefficient is found to be non-zero [19 



In other words, the size of the gluonic contribution is dependent on the definition of the 
parton distributions. In the case where there is a non- vanishing hard gluonic contribution to 
the first moment of gi induced by the anomaly, our analysis would tell us that precisely the 
same contribution will also affect gx- Thus, a measurement of the first moment of grix) will 
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provide a very interesting cross-check about the importance of the anomaly-induced gluonic 
contribution. 

Another interesting prediction for qt could be the analogue of the Bjorken sum rule for 
gi. Given that the first moments of gi and gr are identical, we would expect that gr would 
satisfy a sum-rule which is exactly the same as the Bjorken sum-rule, and whose numerical 
value is the same as that for gi. 

In conclusion, wc have shown that the transverse structure function grix, Q'^){gi{x, Q"^) + 
g2{x,Q'^)) contains a simple operator structure which renders one to understand the spin 
structure of the proton from a completely different experiment involving transversely po- 
larised proton. In addition, due to the simplicity in its structure, it is easier to extract 
and hence understand the higher twist effects. We have shown that the first moment of 
gri^i Q^) measures the spin contributions coming from various partons to the proton spin 
using the Factorisation method. The interesting point to observe is that at large to order 
as{Q'^), with appropriate operator definitions for transverse partons inside the transversely 
polarised proton, the factorisation of mass singularities works. We have found that the 
gluonic contribution to grix, Q^) is zero to this order for the operators discussed. 

Acknowledgements: 

PM and VR would like to thank G. Ridolfi and Diptiman Sen for clarifying some of the 
techniques used in the paper. KS would hke to thank R.G. Roberts for useful discussions. 
VR thanks M.V.N. Murthy for his constant encouragement. We thank the organisers of 
WHEPP-IV for warm hospitahty where part of this work was done. 



13 



References 



[1] D. Adams et al. Phys. Lett. B 336 (1994) 125; K. Abe et al. Phys. Rev. Lett. 
76 (1996) 587. 



[2] 
[3] 
[4] 
[5] 
[6] 
[7] 



R.L. Jaffe, MIT-CTP-2518, |hep-ph 9603^22. 



A.H. Mueller, Phys. Lett. B308 (1993) 355. 

R.P. Feynman, Photon-Hadron Interactions, Benjamin, Reading, 1972. 



R.L. Jaffe, MIT-CTP-2506, |hep-ph 9602236 



R.L. Jaffe, Comm. Nucl. Part. Phys. 14 (1990) 239. 

M. Anselmino, A. Efremov and E. Leader, Phys. Rep. 261 (1995) 1. 



[8] M.J. Alguard et al. Phys. Rev. Lett. 41 (1978) 70; G. Baum et al, Phys. 
Rev. Lett. 51 (1983) 1135; V.V. Hughes et al, Phys. Lett. 212 (1988) 511; 
J. Ashman et al., Nucl. Phys. B328 (1989) 1; B. Adeva et al, Phys. Lett. B 
302 (1993) 533; B. Adeva et al, Phys. Lett. B 320 (1994) 400; D.L. Anthony 
et al, Phys. Rev. Lett. 71 (1993) 759. 

[9] H. Burkhardt and W.N. Cottingham, Ann. Phys. 56 (1970) 453. 

[10] S. Wandzura and F. Wilczek, Phys. Lett. B72 (1977) 195. 

[11] A. Efremov, E. Leader and O. Teryaev, [hep-ph/9607217 . 



[12] R.L. Jaffe, X. Ji, Phys. Rev. D43 (1991) 274. 

[13] J.C Collins, D.E. Soper and G. Sterman, Advanced series on Directions in 
High Energy Physics -Vol 5 (Perturbative QCD) ed. A.H. Mueller. 

[14] G. Aharelli, B. Lampe, P. Nason and G. Ridolfi, Phys. Lett. B334 (1994) 187; 
J. Kodaira, S. Matsuda, T. Uematsu and K. Sasaki, Phys. Lett. B345 (1995) 
527; R. Mertig and W.L. van Neerven, Z. Phys. C 60 (1993) 489; ERRATUM- 
ibid. C65 (1995) 360. 



14 



[15] G.T. Bodwin and J. Qui, Phys. Rev. D41 (1990) 2755. 
[16] R.L. Jaffe and A. Manohar, Nucl. Phys. B 337 (1990) 509. 
[17] A.V. Manohar, Phys. Rev. Lett. 66 (1991) 289. 

[18] X. Ji, Phys. Lett. B289 (1992) 137; A.V. Behtsky, A.V. Efremov and O.V. 
Teryaev, Phys. Atom. Nucl. 58 (1995) 1333, |hep-ph 9512377 . 

[19] C.S. Lam and B. Li, Phys. Rev. D 28 (1982) 683; A. Efremov and O. Teryaev, 
JINR preprint No. E2-88-287 (1988), unpubhshed; G. Altarelh and G.G. Ross, 
Phys. Lett. B 212 (1988) 391; R.D. Carhtz, J.C. Colhns and A.H. Mueller, 
Phys. Lett. B 214 (1988) 229; E. Leader and M. Anselmino, Santa Barbara 
Preprint NSF-88-142 (1988), unpublished. 



15 



